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ABSTRACT 



Aims. The main goal of this study is to perform a sub-arcsecond resolution analysis of the high-mass star formation region 
G 19.61-0.23, both in the continuum and molecular line emission. While the centimeter continuum images will be discussed in 
detail in a forthcoming paper, here we focus on the (sub)mm emission, devoting special attention to the hot molecular core. 
Methods. A set of multi-wavelength continuum and molecular line emission data between 6 cm and 890 yum were taken with the 
Very Large Array (VLA), Nobeyama Millimeter Array (NMA), Owens Valley Radio Observatory (OVRO), and Submillimeter Array 
(SMA). These data were analyzed in conjuction with previously published data. 

Results. Our observations resolve the HMC into three cores whose masses are on the order of 10 1 - 10 M Q . No submm core presents 
detectable free-free emission in the centimeter regime, but they appear to be associated with masers and thermal line emission from 
complex organic molecules. Towards the most massive core, SMA1, the CH3CN (18a- - 17*) lines reveal hints of rotation about the 
axis of a jet/outflow traced by H 2 maser and H 13 CO + (1~0) line emission. Inverse P-Cygni profiles of the 13 CO (3-2) and C ls O (3-2) 
lines seen towards SMA1 indicate that the central high-mass (proto)star(s) is (are) still gaining mass with an accretion rate > 3 10~ 3 
M G yr" 1 . Due to the linear scales and the large values of the accretion rate, we hypothesize that we are observing an accretion flow 
towards a cluster in the making, rather than towards a single massive star. 
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1. Introduction 



individual (G19. 61-0.23) — jets and outflows — Stars: early-type — Submillimeter: 



Understanding the formation of high-mass (M» > 8M Q ) stars, 
as well as their evolution, is the key to investigate the origin of 
the diversity of stars and stellar clusters in the Galaxy. This is 
because high-mass stars not only affect the evolution of the in- 
terstellar medium in general, but are also preferentially born in 
clusters containing a large number of low-mass stars, like, e.g., 
in the Trapezium cluster in Orion (Mu ench et al. 200"2l . Despite 
their importance, the formation process and early evolution of 
OB stars is still far from being understood. In the past decade, 
observations have suggested the existence of flattened rotating 
disk-like structures around newly formed massive (proto)stars. 
In particular, Keplerian disks appear to exist around B-type stars, 
whereas only massive toroids have been found in association 
with O-type stars. These toroids are very different from circum- 
stellar disks, as they appear to be transient, non-equilibrium en- 
tities (Cesaroni et al. 20071 1. In spite of the substantial difference 
between the two types of objects, their existence suggests that 
high-mass star formation may also proceed through accretion 
with angular momentum conservation, as well as in the case of 
low-mass stars. 

This scenario implies the existence of infalling gas onto 
the newly formed stars. Detection of gas infall in the vicinity 
of high-mass stars is hampered by confusion with other pro- 
cesses such as outflow and rotation. Nonetheless, to date inter- 
ferometric observations at centimeter and (sub)millimeter wave- 



lengths have succeeded in detecting infall in a handful of cases 
(e.g., IKeto et al. 19881 ICesaroni et al. 19921 IBeltran et al. 20061 
Zapata et al. 2008 Girar t et al. 20091 1 lending support to the ac- 



* Preprint with full-resolution figures is available at 
http://subarutelescope .org/staff/rsf /publication. html 



cretion scenario. All these findings, however, do not suffice to 
explain the clustered mode of OB-type star formation, which 
is another crucial issue in high-mass star formation theories. A 
way to shed light on this process is to perform detailed observa- 
tional studies of selected star forming regions containing a large 
number of very young OB-type (proto)stars, possibly in different 
evolutionary phases. Excellent signposts of these are ultra com- 
pact (UC) H11 regions and hot molecular cores (HMCs) (see e.g. 
IKurtz et al. 20001 . 

With all this in mind, we performed a detailed study 
of the G 19.61-0.23 high-mass star forming region (SFR) 
which is known to harbor several UC H11 regions and a 
HMC (Garay et al. 1985, 1998; IWobd & Churchwell 19891 
Forster & Caswell 20001 [Furuya et al. 2005a] hereafter paper I). 
The HMC was firstly identified by Garay et al. (1998) in the 
NH3(2,2) inversion transition. A summary of molecular line ob- 
servations towards the HMC up to 2004 is given in Remijan et 
al. (2004; see references therein). No unambiguous proof of the 
presence of OB-type (proto)stars inside the HMC, such as the 
detection of free-free emission, has been found yet. However, 
the det ection of H 2 Q (fHofner & Churchwell 19961 hereafter 
HC96; IForster & Caswell 20001 1 and OH ( |Garayetal. 1985) 
maser emission is considered evidence of their existence. All 
these features make the G 19.61-0.23 star forming region an 
ideal target to study the early phase of the high-mass star for- 
mation process. 
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Table 1. Summary of Interferometric Continuum Emission 
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^major X Utamar 
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Note 


(GHz) 




(U) 


(arcsec) 


(arcsec) 


(deg) 


(mJy beam -1 ) 




4.860 


VLA-A, C 


0.424 - 596.7 


485 


0.630x0.470 


1.9 


0.78 


this work 


8.415 


VLA-A, B, C, D 


0.78- 1027.8 


264 


0.329 x 0.228 


24.4 


18.1 


paper I 


14.940 


VLA-C, B, CnB 


7.72-554.0 


26.7 


0.550x0.500 


-8.4 


1.26 


paper I 


22.272 


VLA-B, D 


1.96-826.8 


105 


0.357 x 0.260 


-4.3 


0.47 


this work 


43.340 


VLA-D 


3.98- 135.6 


51.8 


1.870 x 1.300 


-16.1 


2.57 


this work 


90.700 


OVRO^+NMA^ 


2.90- 139.0 


69.7 


1.52x 1.51 


-9.6 


2.72 


paper I 


335.416 


SMA d 


13.3-241.5 


15.5 


0.85x0.78 


-32 


34.9 


this work 



" The largest detectable size scale. 

6 L, E, H, and UH configurations. 

D, C, and AB configurations. 

d Extended and Compact configurations. 



Based on lower quality images of the region, in paper I, we 
estimated the lifetime ratio between the UC Hn and the HMC 
phases, concluding that the former should last ~3 times longer 
than the latter. With the new observations presented here, we 
wish to set tighter constraints on the statistical study of paper I 
and improve our knowledge of the HMC. In this paper we will 
focus on the latter issue, while a forthcoming paper will be de- 
voted to a more general analysis of the high-mass young stellar 
cluster in the region. 

It is worth noting that the G 19.61-0.23 SFR was believed 
to be located at the near kinematical distance, based on Ho- 
emission and H2CO absorption line observations towards the 
UC Hn regions. This is why in the literature one finds dis- 
tance estimates of 3.8 kpc ( |Georgelin & Georgelin 1976"} , 4.5+1 
kpc ( IDownes et al. 1980b , and 3.5 kpc (Churchwell, Walmsley & 
Cesaroni 1990). However, interferometric observations of the HI 
line at 21 cm seen in absorption against the Hn regions (Kolpak 
et al. 2003) have established that the region is located at the 
far kinematical distance of 12.6+0.3 kpc. A similar result has 
been obtained by other authors (Pandian, Momjian & Goldsmith 
2008), who derive a value of 11.8+0.5 kpc or 12.0±0.4kpc, de- 
pending on the adopted rotation curves. Note that the latter mea- 
surement was done towards the center of the nearby G 19.61 - 
0.13 region, which is offset by 11' from the HMC position. For 
this reason, in this paper, we prefer to adopt d = 12.6 kpc which 
comes from a measurement made towards the center of the re- 
gion of interest for us (G 19.61-0.23). 

Finally, we point out that recently, Wu et al. (2009) have 
reported on the detection of inverse P-Cygni profiles in 13 CO 
7 = 3- 2 and CN 7 = 3- 2 lines towards the G 19.61-0.23 
HMC. This result was obtained using the Submillimeter Array 
(SMA) archive data originally taken by us. Wu et al. (2009) con- 
clude that such profiles are due to infall motions inside the HMC. 
In this paper, we improve on this result presenting a larger study 
of the line emission from the HMC. 



2. Observations and Data Reduction 

Aperture synthesis observations of the continuum and molecu- 
lar line emission towards the G19. 61-0.23 star forming region, 
from centimeter (cm) to sub-millimeter (submm) wavelengths, 
were carried out using the Very Large Array (VLA) of the 



National Radio Astronomy Observatory (NRAOfl (Sect. 12. U . 
the Nobeyama Millimeter Array (NMA) of the Nobeyama Radio 
ObservatoryQ (Sect. 12.2b . the Owens Valley Radio Observatory 
(OVROfl millimeter array (Sect. 12.3b . and the Submillimeter 
An-ay (SMAfl (Sect. |23}, in the period from 2002 to 2007. We 
summarize our continuum imaging and spectral line observa- 
tions in Tables[TJand|2] respectively. 

2.1. VLA Observations 

We performed VLA observations of the continuum emission at 
6 cm, 1 .3 cm, and 7 mm as summarized in Table Q] For all the 
observations described below, we used quasars J1832-105 as a 
phase-calibrator and J1331+305 and/or J0137+331 as flux- and 
bandpass-calibrator(s). 

6 cm — The A- and C-array observations were done on 
November 1, 2004 (project code AF415) and July 11, 2005 
(AF422), respectively, with the standard correlator configura- 
tion for continuum imaging providing a 172 MHz bandwidth 
with dual polarization. 

1 .3 cm — We performed the B-array observations on May 
28 and 30, 2005 (AF415), and the C-array observation on 
November 18, 2006 (AF422). We observed the continuum emis- 
sion in the "BD" intermediate-frequency (IF) pair, with 25 MHz 
bandwidth, and the H2O maser emission in the "AC" pair with 
3. 125 MHz bandwidth and 64-channels, providing a velocity res- 
olution of 0.66 kms -1 . Because of hardware limitations dur- 
ing the EVLA transition phase, we manually supplied the sky- 
frequency of the maser line at the beginning of each switching 



1 The National Radio Astronomy Observatory is a facility of the 
National Science Foundation operated under cooperative agreement by 
Associated Universities, Inc. 

2 The Nobeyama Radio Observatory is a branch of the National 
Astronomical Observatory of Japan, operated by the Ministry of 
Education, Culture, Sports, Science and Technology, Japan. 

3 Research at the Owens Valley Radio Observatory is supported by 
the National Science Foundation through NSF grant AST 02-28955. 

4 The Submillimeter Array is a joint project between the Smithsonian 
Astrophysical Observatory and the Academia Sinica Institute of 
Astronomy and Astrophysics and is funded by the Smithsonian 
Institution and the Academia Sinica. 
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Table 2. Summary of Interferometric Molecular Line Observations 
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H 2 maser 


22235.080 


VLA-B 


0.32x0.23 


-9.6 


0.66 


90 


H 13 CO+ (1-0) 


86754.328 


OVRO 


4.00x3.13 


-21 


0.43 


61 


SiO (2-1) v = 


86846.998 


OVRO+NMA 


3.67x2.40 


-3.0 


0.90 


38 


13 CO (3-2) 


330587.9601 


SMA 


0.94x0.83 


-27 


0.40 


340 


C 18 (3-2) 


329330.5453 


SMA 


0.94x0.85 


-22 


0.80 


338 


CH 3 CNM8 K - 17 K 


331071.594 rf 


SMA 


1.45 x 1.38 


-16 


1.0 


145 



" Effective velocity resolution. 
* Typical RMS noise level per velocity channel. 
c CH» 3 CN 18* - 17* lines as well, see Fig. [TT] 
d ForCH 3 CN£ = 



cycle (12 minutes) instead of using Doppler tracking. We ex- 
cluded all the LHCP data taken at the B-array observations due 
to a 180° phase jump over the correlator band. Furthermore, all 
the data taken with five EVLA antennas during the D array ob- 
servations were excluded due to error on the amplitude calibra- 
tions. We have produced continuum images by merging the B- 
and D-array visibility data (Table [1). For the water maser emis- 
sion, the two data sets were not merged due to variability be- 
tween the B- and D-array observations. We thus prepared a a 3D 
data cube of the maser line emission using only the B-array data, 
which allowed us to attain the best angular resolution. 

7 mm — The 7 mm observations were performed on May 8, 
2007 with the D array, employing the fast-switching technique. 
The adopted switching cycle consisted of a 2.0-minute integra- 
tion on the target and a 50-second integration on the calibrator. 
Before making an image, we excluded all the data taken with the 
nine EVLA antennas due to an unknown error on the amplitude 
calibrations, leaving us with only 17 usable antennas. 

2.2. NMA Observations 

The NMA observations towards G 19.61-0.23 were carried out 
in the period from December 2002 to May 2003 with 3 ar- 
ray configurations (D, C, and AB). We observed the SiO (2- 
1) v = line in the upper-side band (USB) using the FX cor- 
relator with 32 MHz bandwidth centered at the line frequency, 
giving a velocity resolution of 0.108 kms" 1 . For the contin- 
uum and line emission, we used the Ultra Wide Band Correlator 
(UWBC) with a narrow 5 12 MHz bandwidth in each sideband. 
Although this correlator configuration should allow us to ob- 
serve a total bandwidth of 1 GHz with the dual side bands, af- 
ter removing the channels affected by line emission, the effec- 
tive bandwidth usable for continuum emission measurements 
was ~450MHz. We used 3C273 as passband calibrator and 
J1743-038 as phase and gain calibrator. The flux densities of 
J1743-038 were bootstrapped from Uranus, and the uncertainty 
is estimated to be 10%. All the data were calibrated and edited 
using the UVPROC2 and MIRIAD packages. 

2.3. OVRO Observations 

The OVRO observations of G 19.61-0.23 were carried out in 
the period from September 2003 to May 2004 in the three array 
configurations E, H, and UH. We observed the H 13 CO + (1-0) 
and SiO (2-1) v = lines in the lower side-band (LSB). For 



the continuum emission, we simultaneously used the Continuum 
Correlator with an effective bandwidth of 4 GHz and the newly 
installed COBRA correlator with 8 GHz bandwidth. The line 
contamination could not be estimated as these correlators do not 
have spectroscopic capabilities. For molecular lines, we used the 
digital correlator with 15 MHz bandwidth and 60 channels for 
the SiO line, and 7.75 MHz and 62 channels for the H 13 CO + 
line. We used 3C273 and 3C 454.3 as passband calibrators, and 
J1743-038 as phase and gain calibrators. The flux densities of 
NRAO530 and J1743-038 were determined from observations 
of Uranus and Neptune. We estimate the uncertainty of the flux 
calibration to be 10%. The data were calibrated and edited using 
the MMA and MIRIAD packages. 

2.4. Combining NMA and OVRO Data: The SiO (2-1) v = 
Line 

We combined the interferometric SiO (2-1) v = data taken 
with the NMA fSect.12721 and OVRO (Sect.O) interferometers; 
the diameters of the element antennas are 10.0 m and 10.4 m, re- 
spectively. For this purpose, we adopted the observational set- 
ups adopted for both data sets were as similar as possible, with 
two important differences, though: the frequency resolution and 
the integration time for a single visibility. We therefore smoothed 
the NMA visibilities every 8 channels, and averaged them so 
to achieve a 3.35-seconds integration time per visibility. We 
verified that the visibilities taken with the two arrays have a 
fairly nice consistency in flux calibration when we compared 
them over the common range of the projected baseline length. 
Subsequently, we produced a 3D cube using task IMAGR of the 
AIPS package with a robust parameter of +1 to find a compro- 
mise between angular resolution and image-fidelity. 

2.5. SMA Observations and Data Reduction 

Aperture synthesis continuum and line emission observations of 
G 19.61-0.23 at 890//m were carried out with the SMA on May 
12, 2005 in the compact array configuration and on July 8, 2005 
in the extended configuration. The shortest projected baseline 
length, i.e. the shadowing limit, was about 11.9 m. This makes 
our SMA observations insensitive to structures more extended 
than 15'.'5, corresponding to 0.95 pc at a distance of 12.6 kpc. 
The SIS receivers were tuned at a frequency of 335.4158GHz 
to observe the 13 CO (3-2), C ls O (3-2), and CH 3 CN (18-17) 
lines in the LSB, and the HC Io O + (4-3) and CN (3-2) lines in 
the USB. Each side-band covers a 2 GHz bandwidth. The at- 
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E.A,[J200O) 



Fig. 1. Continuum emission maps toward the G 19.61-0.23 high-mass star forming region. All the contour levels are drawn with 

±7<x • 2" where <x is the RMS noise level of the image (TableQ]) and n = 0, 1, 2, The wavelength and FWHM of the synthesized 

beam are indicated in the bottom right corner of each panel. See TableQ]for the values of the synthesized beams. 



tained synthesized beam sizes, effective velocity resolution af- 
ter smoothing, and image sensitivity are summarized in TablesQ] 
and|2] 

We used 3C 454.3 as bandpass calibrator, and J1743-038 and 
J1924-292 as phase and gain calibrators. Flux densities of the 
two calibrators were bootstrapped from observations of Uranus 
and Neptune (J1743-038: 1.64 Jy in May and 1.55 Jy in July; 
J1924-292: 3.16 Jy in May and 3.60 Jy in July), and were sta- 
ble within 12% during the observing period. We estimate a final 
uncertainty on the flux calibration of ~ 20%. The data calibra- 
tion was done using the MIR and MIRIAD packages. Following 
standard calibration procedures, we tentatively subtracted the 
continuum emission from the visibility data, then imaged the 
continuum in each side-band. After verifying that the contin- 
uum emission is confined within a compact region at (Aa, A6) 
= (— 1'/2, +2'.' 2) with respect to the phase tracking center (PTC; 
R.A. = 18* 27'" 38.15 s , Dec. = -11°56'39750 in J2000), we 
subtracted the continuum contribution from the visibility data 
by giving the position offset in task UVLIN in MIRIAD pack- 
age. For the final continuum subtraction process, we identified, 
at least, 33 and 24 lines in the USB and LSB, respectively. 
Subsequently, we constructed continuum emission images with 
natural and uniform visibility weighting functions in conjunction 
with the multi-frequency synthesis technique. For the continuum 
imaging, we did not apply another visibility weighting based on 
the system temperature (r sys ) of each element antenna, because 
the r sys information in the USB for the compact array obser- 
vations were not recorded properly. In contrast, the r sys -based 
visibility weighting was applied to the line data in the LSB, for 
which the T sys was properly recorded. Note that the r sys prob- 
lem affecting the USB was not realized by Wu et al. (2009), who 
thus produced the incorrectly weighted, low-resolution contin- 
uum image shown in their Fig. 1. 

Atmospheric seeing was estimated using the continuum 
emission maps of the (point-like) calibrators, because the appar- 
ent angular diameter of the calibrators is related to the size scale 
characterizing the atmospheric turbulence. The time-averaged 
size scale over each observing track should be comparable to the 



beam-deconvolved diameters of the calibrators. By this means, 
we estimate the average seeing for J 1743-03 8 to be 0715 and 
0795 for the extended and compact configuration observations, 
respectively. We found that the seeing towards J 1924-292 was 
twice worse than that towards J 1743-03 8 in the extended array 
observations. In addition, given the fact that the angular distance 
between J1924-292 and the PTC is rather large (21.8°) we de- 
cided not to use this quasar as phase calibrator. Finally, we re- 
constructed all the images of the continuum and line emission 
data. 

3. Results and Analysis 

3.1. Continuum Emission 

3.1 .1 . Overview of the G 1 9.61 -0.23 Star Forming Region 

Fig-Qjpresents the continuum emission maps, including already 
published 3.8 cm and 3 mm images. As discussed in paper I, the 
cm maps, representing free-free emission, show that the region 
contains a cluster of UC Hn regions, i.e., young massive stars. 
One can see that the overall morphology of the cm emission 
does not change much with wavelength. Most of the apparent 
differences are due to different angular resolution and sensitivity 
to extended structures (see TableQ}. The 7 mm and 3 mm maps, 
which have comparable angular resolutions, also show similar 
structures. Such a similarity implies that these two images con- 
tain both thermal dust emission and optically thin free-free emis- 
sion, as discussed for the 3 mm map in paper I. Note that the 
compact 7 mm emission to the east of UC Hn region A (see pa- 
per I for labeling of the UC Hn regions) is very likely an artifact 
due to inadequate sampling of the visibility plane (Sect. I2.U . 
The SMA 890 fim image looks significantly different from the 
cm- and mm maps, as it shows a bright compact source towards 
the HMC position previously reported. No other 890 pm sources 
than the HMC are detected over the SMA Fo V, above a 3cr upper 
limit of 105 mJy beam , corresponding to a brightness temper- 
ature in the synthesized beam of r s b(3cr) = 0.22 K. In the fol- 
lowing, we describe the continuum results obtained at 890 //m 
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Fig. 2. Continuum emission maps at 890 /im taken with SMA toward G 19.61-0.23 with (a) natural weighting and (b) uniform 
weighting. The solid contours start at the 3<x level and increase in steps of +2cr up to the 50% level for a clarity of the maps. The 
dashed contours start at the -3cr level and decrease in steps of -la. The image noise levels are 54.1 and 34.9 mJy beam for 
panels (a) and (b), respectively. The sensitivity of the uniformly weighted map is better than that of the naturally weighted map. We 
believe that this is due to the fact that visibilities have not been weighted taking into account the system temperature of each antenna 
(see Sect l2.51 l. The synthesized beams are shown in the bottom right corners (l'.'38xl'.'22 with P.A.= -44° for natural weighting; 
0'.'85x0'.'78 with P.A.= -32° for uniform weighting). 



Table 3. Results of Submm Continuum Emission Observations 



Name 


Peak Position (J2000) 




c b 






M£ 


R.A.(hh:mm:ss) Dec.(dd:mm:ss) 


(Jy beairT 1 ) 


(Jy) 


(arcsec) 


fee) 


(M B ) 


SMA1 


18:27:38.069 -11:56:37.30 


1.89 


3.2 


2.4 


0.072 


1300 


SMA2 


18:27:38.021 -11:56:38.75 


0.29 


0.52 


1.5 


0.045 


220 


SMA3 


18:27:37.955 -11:56:35.08 


0.25 


0.15 


0.81 


0.025 


60 



" Peak intensity. The uncertainty may be given by the 3cr level of the image noise level, i.e., 0. 1 Jy beam -1 . 
* Flux density, see Sect. l3.1.2l 

Effective diameter of the submm emission, which is given by 2 \[Ajn where A is the area enclosed by the 5<x level contour enclosing each 
source. Notice that, at SMA2 and SMA3, the 5<x level corresponds to the 61% and 71% levels with respect to their peak intensities, respectively. 
d Effective radius of the submm emission. 

e Mass of the core estimated from the flux S^o/jm in Col. 5, assuming f3 = 1 and a dust temperature of 80 K. 



(Sect. 13.1.21 ) and the high-resolution images at cm wavelengths 
(Sect. [TO). 

3.1.2. 890 yum Continuum Emission 

In paper I, we have demonstrated that the 3 mm continuum emis- 
sion towards the HMC is mostly due to thermal dust emission 
with a small contribution from free-free emission. The latter is 
instead insignificant in the newly obtained SMA image at 890 
/on, which is dominated by the dust continuum emission. Fig. [2] 
shows a close-up view of the SMA 890 /mi images produced 
with natural and uniform visibility weighting functions. It is 
noteworthy that the continuum emission in the natural weighted 
map (Fig. [2^) is elongated to both the south and northwest. The 
uniformly weighted map (Fig. |2j)) clearly resolves the emission 
and shows that the elongation is attributed to two additional 
weak sources, whose peak intensities are slightly above the 7<x 
level. These results indicate that there are (at least) three submm 
sources in the region. Hereafter, we will name these SMA1 (the 
HMC), SMA2 (the core to the south), and SMA3 (the core to the 
north-west), as indicated in Fig. |2]3. 

The most intense 890 //m source, SMA1, is located 0'.'32 
northeast of the peak of the 3 mm "dust" continuum emission 



(see Fig. 2c in paper I and Fig.[3]of this work). This means that 
the 890 fim and 3 mm continuum peaks coincide within the er- 
rors. Therefore, in this study we will identify SMA1 with the 
HMC, whose position is obtained from the 890 fim image (Table 
0. The peak intensity (I v ) of the 890 fim continuum emission 
is 1.89 Jy beam 1 , corresponding to T s \, = 31. OK. The spectral 
index (a) between 3 mm and 890 fim is a3mm-890£im > 2.7, im- 
plying a power-law exponent for the dust emissivity, f3, of > 0.7, 
assuming optically thin dust emission and the Rayleigh-Jeans 
approximation. For this estimate, we integrated the 890 /im flux 
over the three sources SMA1, SMA2, and SMA3 and used the 
continuum flux at 3 mm (147 mJy) obtained from the data of 
paper I, where the free-free continuum contribution from the 
nearby UC Hn region was subtracted by extrapolating the 1 .3 cm 
continuum map to 3 mm under the assumption of optically thin 
emission. We stress that the latter flux is affected by significant 
uncertainties due to the method adopted to subtract the free-free 
contribution (see paper I). Therefore, in the following we prefer 
to adopt /J = 1, which is consistent with the value derived above, 
within the uncertainty, and falls in the range /3 = 1-2 found in 
the literature. 

The 890 fim flux densities (S v ) in Table [3] allow us to esti- 
mate the gas-plus-dust mass traced by thermal emission (Md us t)- 
For this purpose we have assumed a dust temperature (7d us t) of 
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Fig. 3. Comparison of the 6 cm (greyscale plus white contour), 3.8 cm (magenta contour), and 1.3 cm (yellow contour) continuum 
emission maps generated with the common minimum UV distance of 13.3 kA (see Sect. 13.1.31 . Contour intervals for the cm maps 
and the submm map are the same as those in Figs. Q] and [2] respectively, where the lcr noise levels are 0.23, 0.18, and 0.36 mJy 
beatrT 1 and the beams 07433x07321, 07262x07195, and 07344x07251, for the 6.0, 3.8, and 1.3 cm maps, respectively. The largest 
of the three beams is shown in the bottom right. The filled red rectangular and the double red circles indicate the peak positions of 
the 3 mm (paper I) and 890 [im (Sect. [3~TT2b continuum emission, respectively. The filled green circles and light-blue triangles show 
the positions of the OH (Garay et al. 1985) and H2O (Hof ner & Church well 1 996) masers, respectively. Labels A and C refer to the 
UC Hn regions (notation as in paper I and references therein). The filled yellow circle associated with the UC Hn region C indicates 
the peak position of the isolated 7 mm emission seen in Fig.QJ close to the HMC. 



80 K (paper I) and have calculated the values of S v by integrat- 
ing the emission in Fig. [2J3 over the regions inside the 5<x con- 
tour levels of the sources. The value Tjust- 80 K was obtained 
from the fit to the continuum spectrum (see paper I), from the cm 
to the mid-infrared regime, including the newly obtained 7 mm 
and 890 //m fluxes. We estimate that the ambiguity in defining 
the boundary between SMA1 and SMA2 causes a ~ 5% error 
on the value of S v . The resulting masses are given in Table [3] 
Note that these values are approximately inversely proportional 
to the dust temperature and are thus affected by the uncertainty 
on this latter parameter accordingly. Since the assumed temper- 
ature is likely to be correct within a factor 2, we believe that the 
estimated masses are also affected by a similar uncertainty. 



3.1 .3. Centimeter Continuum Emission towards the HMC 

To assess the presence of cm emission towards the HMC, we 
reconstructed the VLA 6 cm, 3.8 cm, and 1.3 cm continuum im- 
ages in Fig.[3]by making use of visibilities whose minimum spa- 
tial frequency range is set to be the same for the 3 bands. We 
chose the shortest baseline equal to 13.3 kA, namely that of the 
890 jt/m maps (see Table [TJ. This allows one to resolve out the 
extended emission. Fig. [3] clearly shows that no free-free emis- 
sion is detected towards the peak positions of the three 890 //m 
sources. Note that SMA1 is located at the center of the H2O and 
OH masers' distributions, whereas no H2O and OH maser spot 
appears to be associated with SMA2 and SMA3. The upper lim- 



its (3cr) obtained from the cm images correspond to brightness 
temperatures over the synthesized beam (T s \,) of 199, 230, and 
30.9 K, respectively at 4.86, 8.42, and 22.27 GHz (see the cap- 
tion of Fig.[3]for the corresponding noise levels and beam sizes). 

We point out that the 7 mm emission seen to the north of the 
HMC (see Fig.QJ) coincides with the cometary UC Hn region C 
(Fig.0. We thus argue that free-free emission from this UC Hn 
region would significantly contribute to the corresponding 7 mm 
continuum emission, whose peak intensity is 25.3 mJy beam -1 
(~10 cr). It is also interesting to note that SMA3 lies right in front 
of the vertex of the cometary shaped UC Hn region C, suggesting 
a physical connection between the two. This could explain the 
cometary shape, with the existence of dense material preventing 
expansion of the ionized gas towards north-west. 

3.2. H 2 masers 
3.2.1. Overall Results 

Fig.|4]shows the H2O maser distribution with respect to the UC 
Hn regions. The color image is a map of the first-order moment 
of the maser lines. We detected two out of the six "maser fea- 
tures" previously reported by HC96, i.e. features N. 1 and 2 in 
their notation. Here, we use the term "maser feature" for a well- 
defined, spatially isolated group of maser spots (see, e.g., HC96 
and |Furuya et a l. 2005b). With an angular resolution of 073, one 
cannot distinguish the maser spots associated with all the differ- 
ent lines seen in the spectrum (Fig. [5). The fact that no maser 
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Fig. 4. First moment map of the H2O masers (color image) and maps of the continuum emission (color contours). The black contours 
overlayed on the maser velocity map correspond to the integrated intensity of the maser spots, with levels starting from the 5cr level 
with 5<t steps. The black, cyan, and magenta contours represent the first 3 contours of the 890 fim, 6 cm, and 3.8 cm continuum 
emission, respectively, (the same as in Fig. [2j. The yellow triangles and green-filled circles show, respectively, the positions of the 
H2O dHofner & Church well 1 996) and OH ( |Garay et al. 1 985) masers. The numbers associated with the yellow triangles are used 
to identify the maser features as done by Hofner & Churchwell (1996). 



emission was detected towards the other four features identified 
by HC96, must be due to the high-variability of H2O masers, be- 
cause the intensities measured by HC96 are all well above our 
3cr sensitivity of 135 mJy beam -1 . 

Figs. |4] and [5] clearly show that feature 1 presents only 
blueshifted emission, whereas feature 2 is mostly redshifted, de- 
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Fig. 5. H2O maser spectra toward the maser features No.l and 
2 obtained by integrating the emission inside the 5<x contour of 
each emission (see Fig.|4|. The vertical dashed line indicates the 
systemic velocity (V sys ) of the cloud, VLSR=41.6km s _1 . 



spite the presence of a few blueshifted lines. Notwithstanding 
the well known high variability of water masers, we note that 
these two features have persisted with approximately the same 
spectral shape since December 1991, when the VLA observa- 
tions of HC96 were made. The terminal velocities (V t ) at which 
blue- and red-shifted maser emissions are seen towards features 
1 and 2 differ by ~ 15-20 km s with respect to the systemic 
velocity (V sys ) of 41.6 km s _1 . All this suggests that the maser 
emission may be originating in a bipolar jet driven by a putative 
young stellar object (YSO) inside the HMC. We will come back 
to this issue in Sect. 13.3.41 



3.2.2. Origin of the H 2 Maser Features 1 and 2 

Features 1 and 2 (see Figs. |4]and|5]l seem to have persisted more 
than a decade. From their distribution (see Fig. 3 of HC96), one 
can argue that only these two features out of the six reported by 
HC96 are associated with SMA1. 

Interferometric observations of H2O masers at 22 GHz 
indicate that they are likely excited in shocked regions at 
the interface between (proto)stellar jets and the ambient 
gas (e.g., |Furuya et al. 2 000), although in some objects these 
masers have been suggested to be tracing rotating disks (e.g., 
ITorrelles et al. 1996b . If the latter were the case of the masers in 
SMA1, from the separation between the blue- and red-shifted 
features (1'.'70) and their relative velocities (~35kms _1 ) one 
could estimate the mass needed to ensure centrifugal equilib- 
rium: this is ~ 3700 M , which is ~3 times greater than that ob- 
tained from the submm continuum emission (see Table [3]). This 
suggests that rotation is not a viable explanation for the kine- 
matics of the masers in the HMC. Moreover, as we will see in 
Sect. 13.3.41 the H 13 CO + (1-0) line appears to trace a bipolar out- 
flow oriented SE-NW. In this scenario the two H2O maser fea- 
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Fig. 6. Total integrated intensity maps of the 13 CO (3-2) (left panels) and C ls O (3-2) (right panels) lines produced with natural 
(upper panels) and uniform (lower panels) visibility weighting functions. Note that the upper panels show the whole area of the 
G 19.61-0.23 star forming region, while the lower panels magnify the region centered on the HMC (greyscale plus thin contour; 
see Table[2]i. For a comparison purpose, the 7<x level contour of the 6 cm continuum (i.e., free-free) emission as in Fig.QJ is shown 
by thin green contours. The magenta stars in the lower-panels indicate the peak positions of the three 890 pm continuum sources 
(Table[3]l. The contours of the CO isotopomer maps start from the 5<x level in steps of 3<x, where the Icr noise levels are 1.9 and 1.8 
Jy beam 1 km s for the natural and uniform weighted 13 CO maps, and 0.8 Jy and 1.7 beam 1 km s for the natural and uniform 
weighted C 18 maps, respectively. We have integrated the emission over the LSR-velocity ranges 30.2 < VLSR/km s < 42.2 for 
13 CO, and 34.7 < VLSR/km s _1 < 44.5 for C ls O. All the other symbols are the same as those in Fig.Q] The emission seen in the top 
of the upper-panels is an artifact due to cleaning problems. 



tures can be associated with a bipolar jet feeding the outflow. We 
thus conclude that the jet interpretation is more likely. 

3.3. 13 CO (3-2), C 18 (3-2), H 13 CO + (1-0) and SiO (2-1) 
Line Emission 

3.3.1. Maps of 13 CO (3-2) and C 18 (3-2) Line Emission 

Fig. [6] presents total integrated intensity maps of the 13 CO and 
C I8 J =3-2 lines made with natural and uniform weightings. 
For the sake of comparison, we plot the 7cr contour level of 
the 6 cm continuum emission map, i.e. the lowest level from 
Fig. ID The natural weighted maps show that the I3 CO and C I8 
(3-2) emitting regions are compact and centered on the HMC, 
while these lines are not detected towards SMA3. Since single- 
dish observations of similar objects indicate that these CO iso- 



topomers trace extended regions, it is reasonable to argue that in 
G 19.61-0.23 the extended emission is filtered out by the inter- 
ferometer. We also note that the 13 CO and C 18 (3-2) lines are 
not seen towards SMA1 in the uniform weighted maps. 



3.3.2. Spectra at the Peak Position of the HMC 

Fig. [7] shows the molecular line spectra obtained from the in- 
terferometric observations towards the peak position of SMA1. 
The H 13 CO + line peaks at V sys , but its line profile is not a sin- 
gle Gaussian. The SiO line shows prominent high velocity wing 
emission. Noticeably, the 13 CO (3-2) and C ls O (3-2) lines do 
not show prominent wing emission, but present deep absorption 
features. In both CO isotopomers redshifted absorption (i.e. an 
inverse P-Cygni profile) is seen, albeit fainter in the C 18 line. 
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Fig. 7. Spectra of four molecular lines in r s b scale towards the 
peak of the 890 /vm continuum emission, i.e., SMA1. The SiO 
and H 13 CO + intensities are scaled by a factor 2.0 in the T s \, scale. 
The vertical dashed line indicates the systemic velocity (V sys ) of 
the cloud, VLSR=41.6km s . In the CO isotopomers spectra, 
we mark the positions of all possible (molecular) lines, calcu- 
lated assuming that they are emitted with velocity equal to V sys 
(see Sect. 13.3- 2b . The thick horizontal blue and red bars denote 
the velocity intervals over which the line emission has been inte- 
grated to produce the maps in Figs. [9] [10] and [16] The emission 
seen around Vlsr - 70 km s _1 of the C 18 spectrum is likely 
due only to the HNCO line. 



The straightforward interpretation of these finding is that the 
HMC is undergoing infall. The LSR-velocity of the deepest ab- 
sorption channel is VLSR=44.5km s _1 for 13 CO, and 46.0 km s _I 
for C ls O, the latter being redshifted by 4.4 km s _1 with respect 
to y sys . The deepest absorption channels of the 13 CO and C ls O 
lines have -29.3 K and -15.0K in T s b, respectively. Note that 
the brightness temperature of the 13 CO line (and even more so 
for the C 18 line) is, in absolute value, less than that of the con- 
tinuum peak (see Sect. 13.1. 2k this indicates that the absorption 
line is real and not an artifact due to resolving out the extended 
emission. 

Finally, we note that we cannot exclude that part of the 
blueshifted wing emission of the 13 CO (3-2) line might be due to 
the CHi, 3 CN K = 7 line, as we will argue in Sect. 13.41 After con- 
sidering all the other molecular lines that could overlap with the 
two CO isotopomer lines, we conclude that line-contamination 
of the 13 CO line over the velocity range between Vlsr ~ 34 
km s" 1 and 50 km s _1 , and of the C 18 line over the range from 
~ km s _1 to 65 km s _1 , is unlikely. 




-20 20 40 _ 1 60 80 100 
V LSR (kms ) 



Fig. 8. Comparison of the absorption features seen in the 13 CO 
(3-2) and HCO + (1-0) lines toward the peak of the 890 yum con- 
tinuum emission. For the sake of comparison, the interferomet- 
ric 13 CO (3-2) spectrum has been obtained after reconstructing 
the data with a synthesized beam equal to the single-dish beam 
of the HCO + line observations with IRAM 30-m telescope (see 
Sect. l3~33T l. 



3.3.3. Comparisons of the Absorption Features Seen in the 
13 CO (3-2) and HCO + (1-0) Lines 

To assess the origin of the redshifted absorption seen in the 13 CO 
and C ls O (3-2) spectra, we made Fig. [8] where we compare the 
I3 CO(3-2) spectrum with an HCO + (1-0) spectrum taken with 
the IRAM 30 m telescope (R. Cesaroni, unpublished data). The 
HCO + spectrum shows a number of prominent absorption fea- 
tures. To make the comparison as consistent as possible, we re- 
constructed the 13 CO image from the visibility data adopting the 
same beam as the HCO + line (#hpbw= 29")- We then took the 
I3 CO spectrum towards the same position observed in the HCO + 
line. Notice that all the absorption features seen in the two lines 
are redshifted with respect to V sys . However, absorption is de- 
tected at different velocities in the two tracers: the absorption 
dips are seen at Vlsr - Vsys = +3.5 km s _1 for 13 CO, and at 
Vlsr - Vsys - +10, +29, and +57 km s" 1 for HCO + . As for the 
two features seen at the highest velocities in HCO + , we believe 
that these are due to clouds along the line-of-sight, because their 
velocities correspond to those of the 21cm HI lines observed 
by Kolpak et al. (2003) - see their Fig. 4. This suggests that all 
the HCO + absorption likely occurs against the UC Hn regions 
in G 19.61-0.23. In contrast, no HCO + absorption is detected 
in the velocity range where 13 CO and C I8 absorption is seen. 
These facts suggest that the latter is due to the HMC, i.e. has a 
local origin, and is thus indicating that the core is undergoing 
infall. The nature of the absorption will be further discussed in 
Sect. 14.2.11 where we study the gas infall in the core. 

3.3.4. H 13 CO + (1-0) Emission 

Here, we present the results of the H I3 CO + (1-0) line observa- 
tions made with the OVRO array (Table[2]l. Fig.[9]presents maps 
of both the bulk and line wing emission. One can see that, while 
the former outlines a structure elongated approximately NE-S W, 
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Fig. 9. (a) Integrated intensity map of the H 13 CO + (1-0) bulk emission. The velocity range used for the integration is 35.5 < 
^LSR/kni s < 48.9. The interval between the thin-black contours is 2<x with the lowest contour corresponding to the 3cr level, 
where the lcr RMS noise is 0.53 Jy beam 1 km s _1 . The thin-white contours are the 95% and 90% levels of the corresponding peak 
intensities, where the 90% levels correspond to the 7 .5<x level. All the other symbols are the same as those in Fig. [6] and the numbers 
associated with the three maser spots to the west are the same used by Hofner & Churchwell (1996). The dashed box indicates the 
area shown in the right-hand panel, (b) Overlay of the blue- and redshifted wing emission maps of the H 13 CO + (1-0) emission. The 
contours are 2<x steps starting from the 3<x level. The blue- and redshifted wing emission maps were obtained by averaging the wing 
emission over the intervals 35.3 < VLsii/km s~'< 37.9 and 43.1 < VLSR/km s~'< 49.2, and their RMS noise levels are 43.8 and 30.2 
mJy beam -1 , respectively. 



the latter can be interpreted as a bipolar outflow oriented SE- 
NW. 

Comparison with the positions of the three submm contin- 
uum sources shows that SMA1 is the closest to the geometrical 
center of the outflow and is hence the most likely candidate for 
powering it. Here the geometrical center is defined as the middle 
point on the line connecting the peaks of the blue- and redshifted 
lobes. We obtain (projected) angular distances to the center of 
3'.'2, for SMA1, and 3'.'8, for SMA2. That SMA1 is at the origin 
of the flow is also supported by comparison with the H2O maser 
jet associated with SMA1, discussed in Sect. 13.21 This jet is also 
oriented parallel to the H 13 CO + outflow and has the same red- 
blue symmetry, which strongly suggests a common origin for the 
two. 

Using the bulk emission map, we have estimated the H 13 CO + 
abundance relative to H2 (X H n co +) from the ratio between the 
mean H 13 CO + column density over the 3cr contour level of the 
890 fj.m continuum emission (see Fig.|9^) and the corresponding 
H2 column density obtained from the submm continuum emis- 
sion. Such an estimate is very sensitive to the temperature of the 
gas and dust. In our calculation, we have assumed that the gas 
and dust are well-coupled, thus having identical excitation and 
dust temperatures. For a fiducial value of 80 K (Sect. [3~L2l >. we 
obtain an abundance of ~ 10~ 10 , but one should keep in mind 
that X H i3 CO + spans a range from ~ 3 10~ n to ~ 6 10~ 10 for T 
varying by a factor 2 with respect to the fiducial value. 

Assuming X H n C0+ = 10~ 10 and T = 80 K (Sect. [3~L2b , 
we have computed the outflow parameters by integrating the 
line emission for the blue and red wings (see caption of Fig. [9] 
for their velocity ranges). With these assumptions, we obtain 
a total outflow mass (M\ be) of 3700 M and a momentum of 
14000 M & km s . Note that, even ignoring the error on X^nco*, 



a factor 2 uncertainty on the gas temperature affects by an 
additional factor ~2 the outflow parameters. The dynamical 
timescale (?d) is estimated from the ratio between the diame- 
ter of the lobes and the difference (in absolute value) between 
the systemic velocity and the terminal wing velocity. We ob- 
tain ~ 8 x 10 4 yr, without correcting for the (unknown) out- 
flow inclination. The mass loss rate (Mfl ow ) and a momentum 
rate (F flow ) are thus 0.05 M yr _1 and 0.17 M G km s _1 yr -1 , re- 
spectively. Values that large are very likely overestimated due to 
the various assumptions made. Nevertheless, they indicate that 
the powering source should be as luminous as ~ 10 s L , accord- 
ing to the empirical relationship derived by Beuther et al. (2002). 
Here, we have implicitly assumed that one is dealing with a sin- 
gle star. If this is correct, then such a star must be in a pre-Hn 
region phase, as no free-free continuum emission has been de- 
tected towards the HMC. We will further discuss this in Sect. 14. II 
and Sect. 1431 

3.3.5. SiO (2-1) v = Emission 

Maps of the SiO Emission — From Fig.[7]one sees that the 
Vt of the SiO emission is blueshifted by ~ 11 km s _1 and red- 
shifted by ~ 15 km s with respect to V sys . The presence of 
HV wing emission strongly suggests the existence of a molec- 
ular outflow driven by a YSO in one of the cores. In order to 
analyze the structure of the SiO emitting gas, we produced maps 
of the bulk emission and HV wing emission in Fig.[l0] In panel 
(b), one sees that the HV wing emission is elongated along the 
east-west direction, with the blue lobe lying to east and the red 
lobe to the west. Such a bipolarity indicates that the SiO HV gas 
very likely traces a bipolar outflow, although the lobes do not 
appear very collimated. The fact that the two lobes are largely 



R. S. Furuya, R. Cesaroni, and H. Shinnaga: Infall, outflow, and rotation in the G19. 61-0.23 hot molecular core 



11 




ia h 27 B1 3a*5 aa 6 c 3/5 
R. A (.12000) 

Fig. 10. (a) Integrated intensity map of the bulk emission of the SiO (2-1) v = line (greyscale plus thin contours) observed with the 
OVRO and NMA interferometers. The green contours have the same meaning as in Fig. [6] The bulk emission has been integrated 
over the LSR-velocity range of 31.9 < VLSR/km s _1 < 48.4, and the RMS noise level is 0.28 Jy beam -1 km s . The yellow stars 
represent the positions of SMA1, SMA2, and SMA3 (see Table[3]l. (b) Overlay of the blue- and redshifted wing emission maps of 
the SiO (2-1) v = emission. The contours are spaced by 2cr and start at the 3<x level. The blue- and redshifted wing emission 
maps have been obtained by averaging the wing emission over the intervals 21.1 < VLsii/km s~'< 31.9 and 48.4 < VlsR/km s~'< 
58.0, and their RMS noise levels are 1 1.6 and 13.9 mJy beam 1 , respectively. The hatched ellipses in the bottom right indicates the 
FWHM of the synthesized beam (Table|2]i. 



overlapping each other suggests that the outflow is seen close to 
pole-on. Interestingly, the velocity structure of the SiO outflow 
is very different from that of the larger scale (~20") 13 CO (2-1) 
outflow mapped by Lopez-Sepulcre et al. (2009) with the IRAM 
30-m telescope. The latter has the redshifted lobe to the NE and 
the blueshifted to the SW. We stress that this is not due to the 
different velocity intervals used to produce the outflow maps, 
because the orientation of the SiO outflow does not change us- 
ing the same velocity ranges adopted by Lopez-Sepulcre et al. 
(narrower than those used by us). We conclude that in all like- 
lihood the larger scale flow is originating from another YSO in 
the region. 

What is the source powering the SiO outflow? We believe 
that it is SMA2, despite the small offset between this and the ge- 
ometrical center of the outflow (see Fig. [T0b). because Fig.fTOh 
shows that the peak of the SiO bulk emission is clearly coinci- 
dent with SMA2 (see Fig. [TOJi). 

Note that this implies that (at least) three bipolar outflows 
are present in the region. Beside the large-scale flow mapped by 
Lopez-Sepulcre et al., we have detected two compact outflows: 
one traced by H 13 CO + (see Sect. 13.3.41) and powered by a YSO 
in SMA1, and another traced by SiO and powered by a YSO em- 
bedded in SMA2. Is the latter as massive as the former? Next, we 
attempt to estimate luminosity of the YSO from the SiO outflow 
parameters. 

Properties of the SiO Outflow driven by SMA2 — To derive 
the mass and kinematical parameters of the outflow, we assume 
that the SiO line is optically thin and the SiO molecule is in LTE. 
We adopt an excitation temperature (r ex ) of 20 K ( Appendix lAl 
and an SiO/H2 abundance ratio of 3x 10~ 9 . The latter is estimated 
from the ratio between the SiO column density (obtained from 
the bulk emission - Fig. [TUh ) and the H2 column density (cal- 



culated from the 890 /jm continuum emission towards SMA 2 - 
Table [TJ. We note that that the outflow masses are likely under- 
estimated because of missing-flux filtered out by our interfero- 
metric observations and the uncertainty in defining a boundary 
velocity (Vb) between the outflowing gas and the quiescent am- 
bient gas. 



We find an outflow mass of 90 M Q and a momentum of 
2100 M kms _I . We also obtain fd - 2 10 4 years, implying an 
outflow mass-loss rate Mfl ow ^ 5 x 10~ 3 M yr _I and a mo- 
mentum rate Fq ovi ^0.06 M kms" l yr 1 . Since outflows are 
believed to be momentum driven (Ca brit & Bertout 19921) . the 
latter value may be taken as an indicator of the outflow strength 
and hence of the mass and luminosity of the YSO powering it. 
Using the relationship between ffl ow and YSO luminosity ob- 
tained by Beuther et al. (2002), we estimate that the YSO pow- 
ering the SiO outflow should be as luminous as ~ 3 x 10 4 L . 
This indicates that the powering source is a high-mass star. 



When deriving the outflow parameters as done here with SiO 
and in Sect. 13.3.41 with H I3 CO + , a word of caution is in or- 
der. The difference between the masses of the H I3 CO + outflow 
driven by SMA1 and the SiO outflow from SMA2 amounts to a 
factor ~40. Such a large number may be due to multiple outflows 
unresolved in our observations as well as to the uncertainties on 
the fractional abundances of the two molecules, which are diffi- 
cult to predict. These caveats cast some doubt on the values de- 
rived in our calculations. Therefore, although the basic conclu- 
sion that the two outflows are associated with high-mass YSOs 
is likely to be correct, the outflow parameters reported here must 
be considered with caution. 
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Fig. 11. Interferometric spectrum of the CH3CN and CH3 CN (18-17) lines towards the peak position of the 890 fim continuum 
source, SMA1. The vertical bars above and below the spectrum indicate the rest-frequencies of the K components of the CH3CN 
and CH3 13 CN 7=18-17 transitions. We also indicate other lines that may be blended with the methyl cyanide lines. 




CH 3 CN and CH 3 "CN 18 K 
0.1 pc 



R.A,(J3000) 

Fig. 12. Integrated intensity maps of the Jk = 18^ - 11 k emission of the CH3CN (green contour) and CH3 13 CN (magenta contour) 
lines. The black dashed contour corresponds to the absorbtions region seen in the 13 CO (3-2) line. The thin black contour represents 
the 7<x level of the uniform weighting 890 fim continuum emission map (Fig.[2j3). All the transitions, except the CH3CN K =0,1, 
and 2 components, are not blended with other lines (see Fig.QT]for the spectrum). The contour intervals are spaced by 2cr and start 
from the 3<x level. In the CH3CN maps, we do not plot contours above the 17cr level to prevent saturation of the maps. The l<x RMS 
noise levels are 6.3, 1.4, 1.3, 0.33, and 0.17 mJy beam 1 km s" 1 for the CH 3 CN (18-17) K = to 2 emission, K = 3, K = 4, and 
CH 3 13 CN (18-17) K = 2, and K = 6, respectively. 



3.4. CHiCN 18* - 17* Emission 
3.4.1. Spectrum 

Fig. QT] shows the spectrum of the CH 3 CN and CH' 3 CN J K = 
18 k - 17k lines observed with the SMA towards the 890 jjm 
continuum peak, SMA1 (Sect. I3.1.2I I. High-excitation K com- 
ponents can be seen in this spectrum. In particular, one can rec- 
ognize all lines up to K =10 for CH3CN, and up to K =6 for 
CH^CN, although overlap with transitions from other species 
(some of which are marked in the figure) cannot be excluded. 
Only the CH 3 CN (18-17) K = 3 and K = 4 and the CH' 3 CN 
(18-17) K = 2 and 6 components do not appear to be affected 
by blending. 

The peak intensities of the K - to 4 lines of the CH3CN 
are comparable each other, suggesting that these CH3CN lines 



are optically thick. One can calculate the optical depth of the 
corresponding CH^CN lines from the ratio between the two iso- 
topomers. The values of r for the 13 C substituted species range 
from 0.1 to 0.5 for K < 5. These imply optical depths of 5 to 
29 for CH3CN, assuming a 12 C/ 13 C abundance ratio of 55 at a 
galactocentric distance of 5.3 kpc (Wils on & Rood 19941 1 . As we 
will discuss in Sect l3.4.3"l opacities that large hinder the usage of 
rotation diagrams to estimate the gas temperature and column 
density. 

3.4.2. Maps: Comparison with the Other Lines 

Fig. Q2] shows maps of some CH3CN and CH 3 I3 CN (18-17) 
lines that are not blended with transitions of other molecules. 
For the sake of comparison, in the same figure we also outline 
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Fig. 13. Rotation diagram obtained from the K — 0, 1, 2, 4, 5, and 
6 components of the CH^CN (18-17) lines shown in Fig.fTTl 

the uniform weighted 890 /mi continuum map (Fig.[2j5) and map 
of the redshifted absorption seen in the 13 CO (3-2) emission, 
described in Sect. 13.3.21 The I3 CO (3-2) map corresponds to the 
velocity channel where the maximum absorption is attained, i.e., 
at Vlsr = 44.0 km s _1 . We remind that no emission is seen to- 
wards SMA1 in the integrated intensity maps of the 13 CO and 
C 18 (3-2) lines (Fig.|6]l. From all this, one obtains the follow- 
ing results: (i) CH3CN and CHPCN line emission is detected to- 
wards SMA1 and SMA2, but not towards SMA3; (ii) the peak of 
the methyl cyanide emission coincides with the peak of SMA1; 
(iii) the redshifted 13 CO absorption is seen towards the center 
of the HMC (i.e., SMA1), but is not detected towards SMA2 
and SMA3; (iv) the emitting region of the low /^-transition, e.g., 
K = 0, 1, and 2, has a size similar to that of the the higher 
/f-lines, while in similar objects (e.g., Beltran et al. 2004) the 
emitting region is smaller for higher excitation transitions. 

3.4.3. Excitation Conditions of CH 3 CN 

Methyl cyanide is known to be an excellent temperature tracer 
and rotation diagrams obtained from the different K components 
are commonly used to estimate the gas temperature and column 
density. However, as argued in Sect l3.4.fl most CH3CN(18-17) 
transitions are optically thick thus making the column density in 
the corresponding level a lower limit. To circumvent this prob- 
lem, we have used the optically thin (t < 0.5 - see Sect J3.4. ft 
lines of CH^CN. Despite heavy blending with other transitions, 
we managed to obtain an estimate of the line parameters by 
fitting groups of K components simultaneously, by fixing their 
separation in frequency to the laboratory values and forcing the 
line widths to be the same. In this way we could successfully 
fit the K — 0, 1,2,4, 5, and 6 components and produce the ro- 
tation diagram shown in Fig. [T3] In our calculation, we have 
assumed a source angular diameter of 2'.'4 from Table [3] and 
a 12 C/ 13 C abundance ratio of 55. The best fit gives a tempera- 
ture of 208 ± 41 K and a source averaged CH3CN column den- 
sity of 2.7 x 10 16 ctrT 2 . While the latter is consistent with the 
values quoted by Wu et al. (2009) and Qin et al. (2010), the 
former is significantly less than their estimates of 552±29 K 
dWu et al. 20091 1 and 609±77 K flQin et al. 201Q) . Such a dis- 
crepancy is due to these authors using the optically thick CH3CN 
lines in their calculations, which leads to a severe overestimate 
of the true rotation temperature. 



(a) (,'HjCN (1H-17) K=7 and LH^CW K=f> 




ia*S? m 3B s S 36=0 
K.A.(J2000) 

(b) CH 3 L'N (1H-17) K=8 and CH 3 12 CW K=B 




K.A.(J2000) 

Fig. 14. Maps of the line peak velocity (color) obtained from si- 
multaneous Gaussian fits to the CH3CN and/or CH^CN J K = 
1%k - 17a: lines (see Sect. 13.4.41 for details). The names of the 
lines above each panel show the transitions that have been ana- 
lyzed. The thin contours are a map of the integrated emission of 
the corresponding lines. The contours are in 2<x steps, starting 
from the 3<x level, the thick contours corresponds to the 5<x level 
of the 890 /itn continuum emission shown in Fig. I2J)). The two 
stars mark the peak positions of SMA1 and SMA2 (Sect. 13.1.21 . 



Since the source angular diameter is comparable to the SMA 
beam and the CH3CN line are optically thick, the line bright- 
ness temperature should be similar to the gas kinetic tempera- 
ture of ~200 K. Instead, only ~20 K are measured with the SMA 
(see Fig.fTTl. implying a beam filling factor of 0.1. Although the 
existence of sub-structures due to dumpiness on angular scales 
smaller than the interferometer beam is very likely, such a filling 
factor appears too small, as we do not reveal important fragmen- 
tation of the HMC in our maps. We thus believe that beam dilu- 
tion may explain only in part the low value of T s b and conclude 
that the optically thick CH3CN lines must be tracing the outer 
regions of the core, where the gas temperature is significantly 
less than the temperature measured in the thinner CH^CN tran- 
sitions, originating from the innermost regions. 

3.4.4. Velocity Structure of the HMC Traced by CH 3 CN 
Emission 

To investigate the velocity field of the innermost part of the 
SMA1 core, in Fig. [14], we plot maps of the CH3CN line centroid 
velocity over the HMC. This velocity field was obtained by fit- 
ting simultaneously multiple A'-components with Gaussian pro- 
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files with separations in frequency fixed to the laboratory values 
(see Pearson & Mueller 1996) and line widths forced to be equal. 
The method employed here is described in e.g., Beltran et al. 
(2005) and Furuya et al. (2008). Considering the line-blending 
and high opacity at low-K lines of CH3CN (Sect. [3ATT >. we de- 
cided to fit the CH 3 CN K = 7 and CHi 3 CN K = 5 lines simul- 
taneously, to obtain the map in Fig. 114a. and the CH3CN K — 8 
and CH 13 CN K — 6 lines, for the map in Fig. [T4b . The values 
of the velocity are displayed only inside the area encompassed 
by the 9cr contour level of the corresponding integrated emission 
map. 

Unlike other cases (Beuther et al. 2005), in our study the 
CH3CN lines appear to trace a clear velocity gradient, increas- 
ing from SW to NE, as observed in similar sources (e.g. Beltran 
et al. 2004). In the G 19.61-0.23 HMC, the inferred velocity 
field maps show a coherent pattern, suggestive of the presence 
of systematic motions of the gas. In addition, the two maps show 
a fairly nice consistency. To further confirm these findings, we 
made velocity-field maps by fitting a single-Gaussian profile to 
the weak, optically thin CH^ 3 CN K — 2 component, which is 
not blended with other transitions. Also in this case, the velocity 
map confirms the existence of a velocity gradient and demon- 
strates that this is not due to opacity effects. 

Note that the velocity gradient is almost perpendicular to the 
axis of the H2O maser jet in Fig. [4] and that of the H 13 CO + 
outflow in Fig. |9j). Therefore, our interpretation of the veloc- 
ity gradient is that the HMC is rotating about the jet/outflow 
axis oriented SE-NW, alike the "toroids" imaged by Beltran et 
al. (2006). 



4. Discussion 

4.1. Evolutionary Stage of the Massive YSOs in the 3 
Submm Sources 

Previous studies (e.g., Codella et al. 1997) have shown that H2O 
masers are preferentially coincident with dense molecular cores 
that do not show continuum emission from ionizing gas, imply- 
ing that such cores may harbor massive YSOs prior to the for- 
mation of an Hn region. This is also consistent with the case of 
low- and intermediate mass stars, where H2O masers are known 
to be associated with the youngest evolutionary phases (Furuya 
et al. 2001; 2006). Given the large masses of the three submm 
cores in G 19.61-0.23 - on the order of 10 1 — 10 3 M Q - one can 
hypothesize that each of them could develop an Hn region; the 
fact that, instead, no free-free emission is detected inside the 
cores (Sect. l3.Ql l. as well as the presence of a H2O masers jet in 
SMA1 (Sect. 13.2b . suggest that the YSOs harbored in the cores 
are massive but in an early evolutionary stage. 

It is possible to constrain the properties of the putative Hn 
regions and corresponding ionizing stars embedded in the cores 
using the upper limits obtained from our observations of the con- 
tinuum emission. For this purpose, we use the method illustrated 
by CTC97 (see also Molinari et al. 2000). 

For the sake of simplicity, we arbitrarily assume that each 
submm source develops a single massive YSO, instead of a clus- 
ter. As explained in CTC97, the peak brightness temperature, 
r s b, of a Stromgren Hn region at a given frequency can be ex- 
pressed as a function of the radius, R^, and Lyman continuum 
photon rate, Nj_, y , of the star. In the calculation, we also adopted 
a source distance of 12.6 kpc, and an electron temperature (T e ) 
of 7200 K, corresponding to the mean for the UC Hn regions in 
G 19.61-0.23 (Garay et al. 1998). For a given T s b, one obtains a 
curve like those plotted in Fig. [15] Note that we have not con- 




Log(R s [pc]) 

Fig. 15. Plot of the brightness temperature in the synthesized 
beam of a Stromgren Hn region, as a function of radius and 
Lyman continuum photon rate (see Sect l4.ll for details). The 
curves correspond to the 3cr upper limits on the free-free con- 
tinuum emission measured towards SMA1 at 4.86 GHz (blue), 
8.42 GHz (green), and 22.273 GHz (red). Here the RMS noise 
levels and the geometrical mean of the beam FWHPs adopted to 
plot are 230 //Jy beam -1 at 5 GHz, 180 yuJy beam -1 at 8 GHz, and 
360 |dy beam -1 at 22 GHz. The labels to the right indicate the 
spectral types of the ZAMS stars with the corresponding value 
of A^Ly, according to Panagia (1973). 



sidered the 7 mm image (Sect. 13.1. 3b . because this has resolu- 
tion and sensitivity about 3-5 times worse than the cm images. 
The three curves correspond to 3<x upper limits obtained from 
the VLA continuum maps at 1.3, 3.8, and 6 cm. The permitted 
values of r s b are those lying under each curve. An additional 
constraint is set by the maximum radius of the Hn region, which 
cannot be larger than the core radius, namely 0.025-0.072 pc, 
depending on the core (TableQ. 

If the putative embedded stars are massive, i.e. earlier than 
approximately B0. 5, they must be also very young, as the corre- 
sponding Hn regions cannot be larger than 6 10 -4 pc =130 AU, 
which means that they are basically quenched. On the other 
hand, we cannot rule out the possibility that the stars are later 
than B0. 5, and in this case the Hn regions could be larger and 
optically thin. We believe, though, that the latter possibility is 
less likely given the above mentioned large masses of the cores 
and the signposts of high-mass star formation associated es- 
pecially with SMA1 (i.e. water masers, high temperature, and 
high-excitation lines), and thus conclude that in all likelihood 
the HMC hides OB type stars in a pre-UC Hn region phase. 

4.2. Velocity Structure of the SMA1 Core: Infall plus Rotation 

The velocity structure seen in CH3CN, a high-density, high- 
temperature tracer, is very different from that obtained from the 
I3 CO and C 18 (3-2) lines (see Fig.[T6i>. A comparison between 
the blue- and red-shifted emission in these lines and the ve- 
locity field of the CH3CN transition is shown in Fig. Q~6] The 
most interesting result is that both the deepest absorption and 
the 890 /im continuum peak roughly coincide with the center of 
symmetry (both in space and velocity) of the CH3CN distribu- 
tion. This configuration is strongly suggestive of the HMC to be 
both rotating about a SE-NW axis and collapsing. The fact that 
no hint of infall, i.e. no redshifted (self)absorption is detected in 
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Fig. 16. The color maps are the same as in Fig. [14] while the 
overlayed solid blue and dashed red contours represent, re- 
spectively, the blueshifted emission and redshifted absorption 
seen in the 13 CO (left panel) and C 18 (right panel). To ob- 
tain these plots, we integrated the line over 3 velocity channels 
centered on the deepest absorption and the most intense emis- 
sion, namely over the ranges: 34.2 <VLSR/km s _1 < 35.8 for the 
blueshifted 13 CO emission, 43.8 <V L SR/km s _1 < 45.0 for the 
redshifted 13 CO absorption, and 36.0 <VLs R /km s _1 < 38.4 for 
the blueshifted C ls O emission, 44.0 < V L SR/km s~'< 46.4 for the 
redshifted C ls O absorption. These velocity ranges are outlined 
by the horizontal blue and red bars in Fig. [7] All the contours are 
in steps of ±2<x and start from the +3cr level. The stars mark the 
peak positions of SMA1 and SMA2. See Sect.|43]for details. 



the CH3CN lines can be explained by these being higher-energy 
transitions than the 13 CO and C 18 (3-2) lines. Very likely, the 
absorption observed in the latter arises in the outer, low-density 
regions of the core where the temperature is significantly less 
than that of the high-density gas traced by CH3CN. In the fol- 
lowing, we discuss the details of this scenario and derive some 
physical parameters characterizing the infalling/rotating gas. 

4.2.1. Infall 

The negative features observed in the 13 CO and C 18 (3-2) spec- 
tra are due to absorption against the innermost, optically thick 
part of the core. Unlike the case of absorption at cm wavelengths, 
where the outer regions of a molecular envelope absorb the free- 
free continuum of an embedded Hn region, here we have de- 
tected absorption against the continuum emitted by the core it- 
self. The problem is thus complicated by the fact that gas and 



dust are mixed, so that both line and continuum photons are 
emitted from any point inside the core. However, in the outer 
region of the core (where absorption occurs) the density is prob- 
ably low enough to decouple the line radiative transfer from that 
in the continuum, whereas in the innermost region (where the 
bright continuum emission comes from) the dust optical depth is 
large enough to absorb all line photons. With this in mind, we 
can simplify the problem assuming a spherical core made out of 
two regions: an outer molecular shell enshrouding an inner, opti- 
cally thick dusty nucleus. The temperature increases outside-in, 
so that the nucleus is hotter than the shell. In this configuration, 
the line brightness temperature measured along the line of sight 
through the center of the core can be written as 

TB = Uv(T ex )-T c ](l-e- T ) (1) 

where T c =s 31 K (Sect. l3~L2l is the measured continuum bright- 
ness temperature, t is the line optical depth, and J v is defined 
as J V (T) - hv/{k[exp(hv/kT) - 1]}, with h and k Planck and 
Boltzmann constants, respectively. We have neglected the con- 
tribution of the cosmological background temperature (Teg - 
2.7 K), as this is very small at the 13 CO and C 18 (3-2) line 
frequencies (hv/k - 15.9 K, hence J V (T BC ) a 0.045 K). 

Equation ([TJ can be written for both the 13 CO and the C ls O 
(3-2) lines, taking into account that the values of T B are respec- 
tively -29 K and -15 K and t( 13 CO) = 6.4 t(C 1s O). From the 
two equations one obtains t( 13 CO) ^4.5 and r ex 7 K. For a 
line width of ~6 km s _1 , these imply a 13 CO column density of 
~ 1.8 x 10 17 cm -2 and an Fb column density Nn 2 — 10 23 cm -2 , 
assuming a I3 CO abundance of 1.5 x 10~ 6 . 

Such a low value of r ex is not consistent with the absorbing 
gas being associated with the HMC and suggests that absorption 
could be due to a cold layer of molecular gas located in the outer 
regions of the cloud. Alternatively, the absorbing gas could cover 
only a fraction of the continuum. In this case, one must introduce 
a filling factor < 1 multiplying J v (T ex ) in Eq. ([TJ and the value 
of 7 K becomes a lower limit. 

Assuming that the absorbing gas is indeed associated with 
the HMC, one can roughly estimate the mass accretion rate for a 
constant density distribution: 

M = 47rm li2 RN ii2 Vi B f (2) 

where we take the infall speed (V; n f 4 km s ) equal to the dif- 
ference between the velocity of the absorption dip and the sys- 
temic velocity (Sect. I3.3.2I I. The radius, R, at which absorption 
occurs is difficult to estimate, but it seems clear that only the 
outermost layers of the core can contribute, because the excita- 
tion temperature derived above (7 K) is much less than the HMC 
temperature obtained from CH' 3 CN (208 K; Sect. 13X31 . Since 
I3 CO is likely thermalized, T ex must be very close to the gas ki- 
netic temperature and this implies that the radius corresponding 
to 7 K must be much greater than that of the CH^CN emitting 
core. Therefore, we can only estimate a lower limit on the mass 
accretion rate assuming a (minimum) value of R equal to the ra- 
dius of the HMC, 0.072 pc (Table©. 

We obtained M > 3 x 10~ 3 M yr _1 , consistent with the find- 
ings of Wu et al. (2009). Note, however, that this is a very con- 
servative lower limit as in all likelihood a temperature of 7 K 
pertains to gas layers located much further than 0.072 pc from 
the HMC center. Indeed, comparison with the large outflow mass 
loss rate (Sect. 13.3.4b suggests that the actual infall rate could be 
much larger than the lower limit quoted above. 

Finally, we note that such a large accretion rate is sufficient to 
quench an Hn region even from a star as luminous as suggested 
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Fig. 17. Position-velocity diagram of the CH 3 CN I83 - 17 3 com- 
ponent toward the SMA1 core. The cut is made along a direc- 
tion with P.A. =37°, i.e. along the velocity gradient identified in 
Fig. [14] The contour levels range from 0.61 to 4.1 in steps of 0.6 
Jy beam 1 . The horizontal and vertical bars in the bottom left in- 
dicate the spatial and velocity resolutions, respectively. The solid 
green curves encompass the region of the plot inside which emis- 
sion is expected if the gas undergoes Keplerian rotation about 
and free-fall infall onto a massive star with mass 83 M . 



by the outflow momentum rate, i.e. 10 5 L Q (see Sect. 13.3.4b . This 
corresponds to a zero-age main-sequence 07 star, which requires 
an infall rate > 2x 10~ 5 M yr -1 to quench the corresponding Hn 
region (see Walmsley 1995), much less than the lower limit de- 
rived by us. It must be noted, though, that the system we are 
considering is not spherically symmetric, so that one cannot ex- 
clude that a hypercompact Hn region would form even under 
such extreme conditions. 

4.2.2. Rotation 

As discussed in Sect. 13.4.41 we believe that the SMA1 core is 
undergoing rotation about the SE-NW axis of the water maser 
jet and H 13 CO + outflow. Is this rotation sufficient to stabilize the 
HMC? We know that infall is occurring on a larger scale than the 
HMC, but it is to be understood if gravitational collapse contin- 
ues inside the HMC or, instead, the infalling material attains cen- 
trifugal equilibrium. The mass that can be supported by rotation 
is Mdyn = R V^ ot /G - 83 M , with G gravitational constant and 
y rot rotational velocity at the outer radius of the HMC, R. The 
fact that the HMC core mass obtained from the sub-mm contin- 
uum emission (~1300 M Q ) is much greater than Md yn argues in 
favor of the core to be undergoing gravitational collapse. 

To check if the idea of a an infalling and rotating HMC is 
consistent with our results, we have made a position-velocity 
diagram of the CH3CN (18-17) K=3 line emission, along the 
direction of the velocity gradient (i.e. approximately NE-SW). 
This is shown in Fig.[17l where we have also overlayed a pattern 
representing the maximum and minimum velocities expected at 
a given position from a rotating and collapsing core. We have 
arbitrarily assumed Keplerian rotation and infall with zero ve- 
locity at infinite distance from the HMC center. The pattern in 
the figure corresponds to a central mass of 83 M Q and an outer 
radius of 10100 AU. 

Although the simple scenario depicted here is far from be- 
ing unique, the pattern is consistent with the line emission, once 
the limited angular resolution is taken into account. This shows 
that one cannot rule out the possibility that the infalling material 



is settling onto a centrifugally supported disk in the innermost 
regions of the HMC. 

4.3. Nature of the SMA1 Core 

What is the stellar content of the HMC? Does this consist of one 
(or a few) massive stars, or is the core hosting a cluster of lower 
mass stars? Our findings do not allow to draw any firm conclu- 
sion, as the molecular gas cannot be investigated with sufficient 
resolution and only a loose upper limit can be set on the HMC 
luminosity. With the advent of the next generation of large tele- 
scopes, it will be possible to overcome these problems, but at 
present we can only make speculations. 

We have seen in Sect. I4.ll that no free-free emission is de- 
tected towards the HMC and that this may imply that no star 
earlier than B0. 5 is present in the core or that the star is still in 
a pre-UCHu region phase. If the scenario proposed in Sect. I4.2I 
is correct, we are dealing with a well defined system, consist- 
ing of a massive core undergoing infall and rotating about a 
water maser jet/H I3 CO + outflow. Such a symmetric structure 
and the large mass accretion rate derived suggest that only few 
high-mass YSOs might be located at the center of this system, 
rather than a cluster with a significant contribution from low- 
mass stars. This hypothesis is supported by the comparison be- 
tween the fragmentation time scale, ff rag , and the free-fall time, 

of the HMC. The former can be estimated as the ratio be- 
tween the core diameter and the line width of a typical HMC 
tracer, e.g. CH 3 CN: f ftag = 0.14 pc/10 kms" 1 a 1.4 X 10 4 yr. 
The latter is calculated for a mean H2 density of ~ 4 x 10 7 cm -3 , 
and is % 6 x 10 3 yr. Because of ff rag ~ 2tg, we argue that 
fragmentation could be partially inhibited during the collapse. 

It is worth pointing out that the previous discussion has lim- 
ited validity, as it is mostly based on qualitative arguments and 
neglects the effect of the magnetic field, which might contribute 
significantly to stabilize the core against gravitational collapse. 
Nonetheless, we believe that our findings lend support to the idea 
that the stellar content of the core could be biased towards very 
young, OB-type stars. 



5. Conclusions and Future Perspectives 

We have performed deep continuum imaging from the centime- 
ter to the sub-millimeter regime of the G 19.61-0.23 high-mass 
star forming region. Our observations have confirmed the exis- 
tence of a large number of UC Hn regions, as well as that of a 
HMC. The latter is resolved into three cores, with one of these 
(here called "the HMC") being ~10 times more massive than the 
others. We have also mapped the region in a number of molecu- 
lar lines at 3 mm and 890 yt/m, most of which appear to trace the 
three cores. Star formation is likely to occur not only in the HMC 
(SMA1), but also in the other two (SMA2, SMA3), as witnessed 
by the existence of a bipolar outflow seen in the SiO (2-1) line 
towards SMA2. 

The CH3CN (18-17) line emission reveals a velocity gra- 
dient across SMA1, roughly perpendicular to a water maser jet 
and bipolar H I3 CO + outflow directed in the SE-NW direction. 
We interpret this velocity gradient as due to rotation about the 
jet/outflow axis. We also confirm the existence of an inverse P- 
Cygni profile in the 13 CO (3-2) line, already detected by Wu et 
al (2009) and reveal a similar profile also in the C ls O (3-2) line. 
This redshifted absorption strongly suggests that the SMA1 core 
is infalling - beside being rotating - with a mass accretion rate 
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> 3 10~ 3 M Q yr _1 . We conclude that very young OB-type stars 
are likely forming inside the HMC. 

The study presented here is an excellent benchmark of what 
will be feasible, with much better resolution and sensitivity, 
with new generation instruments. Deep unbiased surveys of se- 
lected high-mass star forming regions will permit to improve our 
knowledge of the OB star formation process, identifying newly 
formed stars by means of their free-free continuum emission and 
investigating at the same time the structure and velocity field of 
the molecular cores associated with them. 
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Appendix A: Rotation Diagram Analysis of the SiO 
Lines 

We estimated the excitation temperature of the SiO molecule 
from the (2-1), (3-2), and (5-4) transitions (Fig. IA~T1 J. M. 
Acord, personal communication) that have been observed simul- 
taneously with IRAM 30 m telescope. In this way, problems due 
to relative gain calibration and pointing errors are minimized. 
Following Acord et al. (1997), we have applied the "rotation di- 
agram" method, whose advantages and limitations are discussed 
in detailed in (Goldsmit h & Langer 1999| l. For this purpose, in 
Fig. IA.2l we plot the logarithm of the column density in the lower 
level of each transition, divided by the corresponding statistical 
weight, against the energy of the upper level, E u . The column 
densities were obtained from the line intensities under the as- 
sumption of optically thin emission. The column densities were 
also corrected for the different beam filling factors of the three 
lines, assuming the source to be point-like and referring all mea- 
surements to the 28" beam at the frequency of the (2-1) line. 

The slope of the linear fit to the data, l/T, gives the "rota- 
tional temperature" (r rot ) of the SiO molecules, which is likely 
an underestimate of the kinetic temperature of the H2 gas (see 
e.g. Acord et al. 1997). The data are well fitted by a straight line, 
which indicates that the method used is likely correct, although 
the SiO molecule is known to be subthermally excited as noted 
above. We obtain r rot -20K and a total SiO beam averaged col- 
umn density of 4.6 10° ctrr 2 . 
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Fig.A.l. Spectra of thermal (v = 0) emission of the SiO 
molecule towards G 19.61-0.23 taken with the IRAM30m 
telescope (J.M. Acord personal communication) in main-beam 
brightness temperature (r m b) scale. The vertical dashed-line in- 
dicates systemic velocity (V S ys) of Vlsr = 41.6 km s _1 . 



29 



26 



1 1 1 1 1 1 1 1 1 1 




i - 


I_ G19.61 


-0.23 




SiO (J+1 


J) v=0 




~_ T M = 20 K 






_ N, ,(SiO) = 4.6E+13crrT 2 






- I I I I I I I I I I I 




i - 



10 20 

E JK A (k) 



Fig. A.2. Rotation diagram of the SiO thermal emission shown 
in Fig. IA. 1 1 The dashed line is the best-fit to the data. The values 
of the rotation temperature and column density thus obtained are 
given in the bottom left of the figure. 



